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^N ' Abstract. A 7.5 years spectrophotometric monitoring program of 28 

Qh' Palomar-Green quasars to determine the size of their broad emission hne 

.^ I region (BLR) is reviewed. We find both the continuum and the emission 

hne fluxes of ah quasars to vary during this period. Seventeen objects 

0^ \ has adequate samphng for reverberation mapping and in aU of them we 

^^ ■ find the Balmer hne variations to lag those of the continuum by ~100 

I \ days. This study increases the available luminosity range for studying 

^ ■ the size-mass-luminosity relations in active galactic nuclei (AGNs) by two 

OO \ orders of magnitude and doubles the number of objects suitable for such 

22 ' studies. Combining our results with data available for Seyfert 1 galaxies, 

^T \ we find the BLR size to scale with the rest-frame 5100 A luminosity as 

Q ■ ^o.70±o.03^ This result is different from previous studies, and suggests that 

^^ \ the effective ionization parameter in AGNs may be a decreasing function 

O ■ of luminosity. We are also able to constrain, subject to the assumption 

r"| . that gravity dominates the motions of the BLR gas, the scaling relation 

Qh' between the mass of the central black hole and the AGN's luminosity. 

We find that the central mass scales with the 5100 A luminosity as M oc 
^ . ^0.5±0.1_ 

(^ . A program to monitor 11 high-luminosity quasars is presented here 

for the first time. Preliminary results from this program indicate contin- 
uum variation of order of 0.1 mag in all objects. We illustrate the impor- 
tance and feasibility of monitoring those objects spectrophotometrically. 
When this program will be completed reverberation mapping studies will 



C^ ' cover the entire AGNs luminosity range from 10^^ to 10^^ ergs s ^. 



1. Introduction 

Broad emission lines in active galactic nuclei (AGNs) emerge from the innermost 
regions of these objects. As such, they provide important information (e.g., 
composition, dynamics, physical conditions, and geometry) about the AGNs' 
unresolved regions. Reverberation mapping, observing the degree and nature of 
the correlation between continuum and emission-line flux variations, is one of 
the major tools for studying the distribution and kinematics of the gas in the 
broad line region (BLR) and to study the central masses of AGNs (e.g., Peterson 
1993; Netzer &: Peterson 1997). At a first approximation reverberation mapping 
yields a measure for the size of the BLR. Combining this size with the line profile, 
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which represent the kinematics, one can estimate the mass of the central source 
in an AGN. Determining this mass is most important for the understanding 
and modehng of the AGNs phenomenon. Estabhshing the relations between the 
physical properties (such as BLR size, mass of the central source, and the AGN 
luminosity) of an AGNs sample will provide us with important information for 
understanding the characteristics which are common to all AGNs. 

In the past ~ 15 years about 17 low-luminosity AGNs (Seyfert 1 galaxies) 
have been successfully monitored and produce statistically meaningful BLR sizes 
(see Wandel, Peterson, & Malkan 1999, and references therein). Best studied 
among these is the Seyfert 1 galaxy, NGC 5548, which was monitored from the 
ground for over eight years, and from space for several long periods (Peterson 
et al. 1999, and references therein). Several other Seyfert Is were observed for 
periods of order 1 year or less, and nine Seyfert Is were studied over a period 
of eight years (Peterson et al. 1998a). The measured time lags between the 
emission lines and the continuum light curves in these objects can be interpreted 
in terms of the delayed response of a spatially-extended BLR to a variable, 
compact source of ionizing radiation. While the observations do not uniquely 
determine the geometry of the BLR, they give its typical size which, for Seyfert 1 
galaxies, is of the order of light-days to several light-weeks (~ 10^^-10^^ cm). 
Recent studies have shown that the time lags determined in NGC 5548 for 
different observing seasons correlate with the seasonal luminosity of the object 
(Peterson et al. 1999), and have presented evidence for Keplerian motions of 
the BLR gas (Peterson & Wandel 1999). 

While there has been great progress in mapping Seyfert I's few similar 
studies of the more luminous AGNs - the quasars - have been presented. There 
still have been some open questions such as, do quasar emission lines respond 
to the continuum changes, as seen in Seyfert galaxies? What is the relative 
amplitude of the response? What is the lag of the response, reflecting the light- 
travel time across the BLR? Do quasar BLR sizes scale with AGNs luminosity, 
and lie on a continuous relation from the faintest Seyferts to the bright quasars? 

There are several difficulties when attempting to monitor high-luminosity 
AGNs. Quasars have fainter apparent magnitudes, hence one needs larger tele- 
scopes and/or much longer integration limes. Since quasars variability time 
scales are longer than Seyferts I's and their BLRs are expected to be an order of 
magnitude larger than in Seyfert I's, we need much longer monitoring periods. 
When monitoring Seyfert I's one often use the narrow emission lines to intercal- 
ibrate the observed spectra, however, in quasars narrow emission lines are very 
faint (or not present at all) and a different relative flux calibration method needs 
to be exploit. 

Past attempts to spectrophotometrically monitor quasars have generally 
suffered from temporal sampling and/or flux calibrations that are not sufficient 
for the determination of the BLR size (e.g., Zheng 1988; Perez, Penston &: 
Moles 1989; Korista 1991; O'Brien & Gondhalekar 1991; Jackson et al. 1992; 
Koratkar et al. 1998; Wisotzki et al. 1998). The quasar best studied by lUE, 
3C 273, has yielded disputed results when different researchers have analyzed 
similar lUE monitoring data sets. Both O'Brien & Harries (1991) and Koratkar 
& Gaskell (1991a) found a measurable and similar lag between continuum and 
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BLR variations, while Ulrich, Courvoisier, & Wamsteker (1993) argued that the 
hne variations reported in the earher studies were only marginally significant. 

As a more definite results on the BLR size in quasars is needed we began two 
quasars' reverberation mapping programs on which we report in this contribu- 
tion. The first program is the monitoring of a sub-sample of 28 Palomar-Green 
(PG) quasars. The program and its results are described in details by Kaspi 
et al. (2000) and references therein. Here we present the basic concepts and 
results of this program and discuss the size-mass-luminosity relations in AGNs. 
The second program of monitoring 11 high-luminosity high-redshift quasars is 
reported here for the first time and some preliminary results are presented. 



2. The PG quasars — Sample & Observations 

In contrast to Seyfert 1 galaxies, which were chosen for reverberation mapping 
mainly because of their variability properties, the quasars in our sample were 
selected according to their observable properties. The sample of 28 quasars is 
drawn from the 114 PG quasars (Schmidt & green 1983). The PG quasars are 
well studied objects over the whole electromagnetic spectrum from X-ray to 
Radio (e.g., Neugebauer et al. 1987; Kellermann et al. 1989; Boroson & Green 
1992; Brandt, Laor, &: Wills 2000). Many of their observational properties are 
well known and adding the information about their variability properties and 
BLR sizes will increase our understanding of the AGNs phenomenon. 

We selected objects with northern declination, B < 16 mag, redshift z < 
0.4 (so that the Balmer lines can be observed in the optical region), and a 
bright comparison star within 3'5 of the quasar. The absolute magnitude range 
covered by the sample is —23 > Mb > —27 mag (using qq = 0.5, Hq = 75 
kms~^Mpc^^, and zero cosmological constant) and the bolometric luminosity 
range is 4 x 10^'^ < L < 3 x 10^^ ergss"^. 

Spectrophotometric observations of the sample were carried out form 1991 
March, over a period of 7.5 years, until September 1999. We used the Wise 
Observatory (WO) 1 m telescope to observe the sample once a month (when the 
objects are observable) and the Steward Observatory (SO) 2.3 m telescope once 
every ~ 4 months. The spectral range is from 4000 to 8000 A with spectral 
resolution of ~ 10 A. The Spectrophotometric calibration for each quasar was 
accomplished by rotating the spectrograph slit to the appropriate position angle 
so that the nearby bright star was observed simultaneously with the quasar. 
A wide slit was used to minimize the effects of atmospheric dispersion. This 
technique provides excellent calibration even during poor weather conditions, 
and accuracies of order l%-2% can easily be archived. 

Alongside the spectrophotometric monitoring we carried out at the WO 1 m 
telescope a broad band photometric monitoring in B and R. The 28 quasars 
in our sample were included in a sample of 45 PG quasars which were moni- 
tored monthly to find the variability properties of the PG quasars (Giveon et al. 
1999). We used the photometric data (by means of differential photometry with 
other stars in each field) to check on the continuum behavior found in the the 
spectrophotometric monitoring and to verify that non of our comparison stars 
are variables to within ~2%. We also used the photometric observations to add 
additional epochs to the continuum light curves of the quasars. 
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Continuum light curves were extracted at rest wavelength of about 5100 A 
and line light curves of Ha, H/3, and H7 (where available) for all objects. The 
photometric data were combined into the continuum light curves. Light curves 
for two quasars are shown in Figure |^. PG 0804+761 is our best sampled object 
with ~ 70 spectroscopic observations and ~ 40 photometric observations. The 
light curves of this object clearly show the variations of the three Balmer lines to 
lag the continuum variations. Also they demonstrate how small variations in the 
continuum light curve smear out in the line light curves - this effect is a result 
of a stratified BLR. PG 1229+204 is a more typical example of our sample. It 
has ~ 30 spectroscopic observations and ~ 30 photometric observations. In this 
object it is hard to identify the time lag by looking at the light curves but the 
use of cross-correlation is clearly yielding a time lag (see next section). 



3. Observational Results and Time Lag Determination 

At the end of the 7.5 years monitoring period we find that for 17 objects out 
of the 28 there are more than 20 spectroscopic observations. Typically there 
are between 20 to 70 observations for each of these 17 quasars. The other 11 
quasars have less than 10 spectroscopic observations (the typical number is 5 
observations) which is not adequate sampling for time series analysis. These 11 
objects are excluded from further discussion hereafter. 

All the 17 quasars with adequate spectroscopic sampling had gone contin- 
uum variation which, quantified as -Fmax/-^min — 1, lie in the range of 35%-150% 
for different objects. All 17 quasars also show line flux variation which follow the 
continuum variation with an amplitude of about half of the continuum variations 
(see Figure ||). 

In order to determine the time lags of the line light curves relative to the 
continuum light curves we use two cross correlation methods: one is the in- 
terpolated cross-correlation function (ICCF; Gaskell &: Sparke 1986; Gaskell & 
Peterson 1987; White & Peterson 1994) where one light curve is interpolated and 
then it is cross-correlated with the second observed light curve; then the second 
light curve is interpolate and cross-correlated with the first observed light curve. 
The final ICCF is the mean of these two cross-correlation functions. The second 
method is the z-transformed discrete correlation function (ZDCF) of Alexan- 
der (1997) which is an improvement on the discrete correlation function (DCF) 
of Edelson & Krolik (1988). This method applies Fisher's z transformation 
to the correlation coefficients, and uses equal population bins rather than the 
equal time bins used in the DCF. The two independent methods are in excellent 
agreement for our data and in the following analyses only the ICCF results are 
used. Figure ^ demonstrate the cross correlations of two of the line light curves 
presented in Figure ffl with their corresponding continuum. For the purpose of 
this work we used the centroid of the ICCF (computed from all points within 
80% of the ICCF peak value) to define the time lag (Gaskell 1994 and references 
therein). 

To determine the uncertainties in the cross-correlation time lag we used the 
model-independent FR/RSS Monte Carlo method of Peterson et al. (1998b). 
In this method, each Monte Carlo simulation is composed of two parts: The 
first is a "random subset selection" (RSS) procedure which consists of randomly 
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Figure 1. Light curves for two PG quasars. Circles are spectrophoto- 
metric data from WO, squares are spectrophotometric data from SO, 
triangles are photometric data from WO. Continuum flux densities, 
fA, are given in units of 10"^^ ergss~^ cm~^ A"-*^. Emission-line fluxes 
are displayed in units of lO^^^ergss"^ cm~^. Horizontal axis given in 
Julian Day (bottom) and UT date (top). 
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Figure 2. ICCF (solid line) and ZDCF (circles with error bars) for 
two cases. The two cross-correlation methods are in excehent agree- 
ment. Top panel: PG 0804+761 H/3 cross correlated with the optical 
continuum; the ICCF yields time lag of At(centroid) = 151^24- Bottom 
panel: PG 1229+204 Ha cross correlated with the optical continuum; 
the ICCF time lag is Af(centroid) = 71^11- 
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drawing, with replacement, from a light curve of N points a new sample of 
N points. After the N points are selected, the redundant selections are re- 
moved from the sample such that the temporal order of the remaining points 
is preserved. This procedure accounts for the effect that individual data points 
have on the cross-correlation. The second part is "flux randomization" (FR) 
in which the observed fluxes are altered by random Gaussian deviates scaled 
to the uncertainty ascribed to each point. This procedure simulates the effect 
of measurement uncertainties. Applying the above procedure to the line and 
continuum light curves and cross-correlating them is considered one realization 
of the Monte Carlo simulations. Using such ~ 10000 realizations builds up a 
cross-correlation peak distribution (CCPD; Maoz & Netzer 1989). The range 
of uncertainty contains 68% of the Monte Carlo realizations in the CCPD and 
thus would correspond to la uncertainty for a normal distribution. Peterson et 
al. (1998b) demonstrate that under a wide variety of realistic conditions, the 
combined FR/RSS procedure yields conservative uncertainties compared to the 
true situation. 

From the 46 line light curves 40 result with a significant correlation (peak 
correlation coefficients above 0.4). All 40 CCFs indicate a positive time lag of 
the Balmer lines with respect to the optical continuum. The time lags are of 
order of a few weeks to a few months, and the CCF peaks are highly significant 
for most lines. We conclude that a time lag has been detected in one or more of 
the Balmer lines for all 17 quasars. 



4. Size, Mass, and Luminosity — Determination and Relations 

To construct the largest sample with available reverberation mapping data we 
analyze our data together with comparable published data for other AGNs. 
Wandel et al. (1999) have uniformly analyze reverberation mapping data of 
17 Seyfert I's and deduce time lags using the same techniques described in the 
previous section. Combining their results with ours we obtain reliable size- 
mass-luminosity relations for 34 AGNs spanning over 4 orders of magnitude in 
luminosity. 

BLR Size: Since we have both Ha and H/3 time lags for many objects we 
average the two to get a better estimate for the Balmer lines time lag. We do 
not include H7 in the mean since its light curves are very noisy due to the small 
S/N of the line, and the uncertainty of the H7 time lag is consistent with zero in 
several cases, hence, counting it in the mean will add noise into our results. The 
BLR size is then computed as the mean of the Ha and H/3 time lags divided by 
a factor of 1 -|- z to account for the cosmological time dilution. For the Seyfert 
I's we use the H/3 time lags from Wandel et al. (1999) and correct them by the 
(1 -|- z)^^ factor. 

AGN Luminosity: A major limitation in the luminosity determination 
is the lack of knowledge about the ionizing continuum and the spectral energy 
distribution (SED) of the objects in question. Much of the ionizing continuum 
is emitted in the unobservable far-UV and there are still unsolved fundamental 
issues concerning the shape of the continuum (e.g., Zheng et al. 1997; Laor 
et al. 1997). Even in one of the best studied AGNs, NGC 5548, the SED 
is poorly determined (Dumont, Collin-Souffrin, & Nazarova 1998). Another 
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complication is the contribution of the host galaxy to the luminosity of the 
nucleus. Since resolving these complications is an issue for an in-depth study 
we took the simplified approach (following Wandel et al. 1999) of using the 
monochromatic luminosity, XLx, at 5100 A (rest wavelength) as our luminosity 
measure (assuming deceleration parameter qq = 0.5, Hubble constant of Hq = 
75 kms~^, and zero cosmological constant). The uncertainty in this value is 
taken to be only due to the variation of each object and is represented by the 
root mean square (rms) of the continuum light curve. 

Central Mass: Estimation of the AGN's central mass is carried out by 
assuming gravitationally dominated motions of the BLR clouds: M ~ G~^v^r 
(e.g., Gaskell 1988; Wandel et al. 1999; Peterson & Wandel 1999). In this 
relation the radius, r, is the BLR size computed above and the velocity, f , is 
estimated from the rest-frame FWHM of the emission line. Because the broad 
emission lines of AGNs are composed of a narrow component superposed on a 
broader components, a unique FWHM determination is not straightforward. We 
took two approaches to measure the FWHM: the first approach is to measure the 
FWHM of the Balmer lines in each spectrum for a given object and then to use 
the mean FWHM, wfwhm (mean) . The second approach is the one proposed by 
Peterson et al. (1998a) of using the rms spectrum to compute the FWHM of the 
lines, VFWHM(rms). In principle, constant features in the mean spectrum (such as 
narrow forbidden emission lines, narrow components of the permitted emission 
lines, galactic absorptions, and constant continuum and broad-line features) are 
excluded in this method. The FWHM from the rms spectrum measures only 
the part of the line that varies and thus corresponds to the BLR size measured 
from the reverberation mapping. In the following we uses the two approaches 
together in order to compare them. 

Following the approach of averaging the Ha and H/3 time lags, we also 
average the FWHM of the Ha and H/3 lines. To calculate the mass we also 
introduced a factor of v3/2, to account for velocities in three dimensions and 
for using half of the FWHM. The virial "reverberation" mass is then: 



,5 f -Rblr \ f ■wfwhm 



2 



M = 1.464 X 10^ f^ ": , Mq . (1) 

Vltdays/ VlO^kms^V ^ ^ ^ 

In Table |l] we present the above computed properties for all 34 AGNs in 
our combined sample. 

4.1. Size Luminosity Relation 

The BLR size versus the luminosity is plotted in Figure y. The correlation 
coefficient is 0.827, and its significance level is 1.7xl0~^. A linear fit to the 
points gives 



/AL.(5100A)V'°°^°-°'' 



i^BLR = (32.91- j (^^^5^^^ j It days (2) 

(solid line plotted in Figure ^). Considering the Seyfert nuclei (log(ALA(5100 
A));^ 44.2) or the PG quasars alone, we find only marginally significant corre- 
lations, probably because of the narrow luminosity ranges. A significant corre- 
lation emerges only when using the whole luminosity range. 
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Table 1. Sizes, Luminosities, & Masses 



Object 


-Rblr 


ALa(5100A) 


M(mean) 


M(rms) 




(It-days) 


10^4 ergs s~i 


IO^Mq 


IO^Mq 
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42t^i 
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2.3+1-? 


3.0+li 
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37115 
18.711-0 
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16.3+?:^ 

1 4+3-3 
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'-'•'' -0.29 
n. 00+0. 14 
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-i q+2-l 
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3.0lli 
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i-2olS:i 
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Figure 3. BLR size-luminosity relation. PG quasars are denoted by 
circles and Seyfert I's are denoted by squares. The solid line is the best 
fit to the data. The dashed line is a fit with a slope of 0.5 . 



The relation we find between the BLR size and the luminosity does not 
agree with earlier studies which found a smaller power-law index (closer to 0.5, 
e.g., Koratkar & Gaskell 1991b; Wandel et al. 1999). A line with 0.5 slope was 
fitted to the data and is shown as a dashed line in Figure ^. The combined 
sample is clearly inconsistent with this slope. 

Also, under the assumptions that the shape of the ionizing continuum in 
AGN is independent of the luminosity, and that all AGNs are characterized by 
the same ionization parameter and BLR density (as suggested by the similar 
line ratios in low- and high-luminosity sources), one expects Rblr oc L^'^. This 
theoretical prediction is based on the assumption that the gas distribution, and 
hence the mean BLR size, scales with the strength of the radiation field. Our 
present result suggests that those assumptions should be re-examined. In par- 
ticular if we keep the assumption that the BLR density is the same for all AGNs 
then the ionization parameter, U, should have the relation U ex L~^'^. 
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4.2. Mass Luminosity Relation 

The mass-luminosity relation is plotted in Figure |^ for the two mass estimates 
described above. Our mass estimates based on the determination of the FWHM 
measured from the rms spectra are plotted in the top panel. The correlation 
coefficient is 0.473 with a significance level of 4.7x10^'^. A linear fit to this 
relation gives 

M(,n.) = (5.75«:i) X 10' (^^^) Afe (3) 

and is plotted as a solid line in the diagram. 

The mass estimates based on the determination of the FWHM from the 
mean spectra are plotted versus luminosity in the bottom panel of Figure 0. 
The correlation coefficient between these two parameters is 0.646 and has a 
significance level of 3.7x10"^. A linear fit gives 

h0.46A ^in7M^A(5100A)^ 



M(mean) = (5.71^--j x 10^ [^J^^^^, j M^ (4) 

and is also plotted as a solid line. 

A surprising result is that when using the FWHM from the rms spectra the 
mass-luminosity relation is less significant than when using the mean FWHM. 
This is in contradiction to the theoretical considerations leading to the use of the 
FWHM from the rms spectra (see above). This disagreement can be attributed 
perhaps to the fact that the line fluxes in the rms spectra are weaker and hence 
the uncertainty in the corresponding FWHM might be larger. 

Weighting the two mass-luminosity relations according to their significance 
our results imply M oc l0.5±o.i^ This does not agree with previous results - 
for example Koratkar & Gaskell (1991b) found M oc L°-^'^^o.25 g^^^ Wandel 
et al. (1999) reported on M oc L°-'^'^'^^-07 _ The fact that the scatter in the 
mass-luminosity relation is larger than that of the size-luminosity relation may 
indicate that luminosity, rather than mass, is the variable that mainly determines 
the BLR size. 

Using a rough estimate for the bolometric luminosity as -Lboi^ 9AL;^(5100A), 
we obtain an Eddington ratio of 

Lboi ^^,.^ ALa(5100A) V-' 

The Eddington limit, based on this rough estimate for Lbob is plotted as a dashed 
line in Figure ^ This result indicates, for the first time from reverberation map- 
ping studies, that the Eddington ratio increases with luminosity. Our findings 
are inconsistent with theoretical models proposing the AGNs' luminosity to be a 
set fraction of the Eddington luminosity, e.g., geometrically thin, optically thick 
accretion disk model which implies M oc L (Laor & Netzer 1989). Our result 
suggests that the mass accretion rate grows with luminosity much faster than 
the central mass, which would mean very different disk properties in low- and 
high-luminosity sources. 
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Figure 4. Mass-luminosity relations. Top: masses derived from 
iTWHM (rms) . Bottom: masses derived from ijfwhm (mean) . Symbols 
as in Figure ^. Solid lines are the best fit to the data. Dashed lines 
are the Eddington limit based on a rough estimate for the bolometric 
luminosity (see text). 
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5. Reverberation Mapping of High-Luminosity Quasars 

Though the PG quasars are two orders of magnitude more luminous than Seyfert 
1 galaxies they are considered to be low-luminosity objects among the known 
quasars. We can enlarge the luminosity range of AGNs with reverberation map- 
ping data by two orders of magnitude if we study high-luminosity, high-redshift 
quasars {L > 10^^ ergs s~^). In fact, nowadays knowledge regarding line vari- 
ability properties of high-luminosity quasars resembles what was known for the 
PG quasars and low-luminosity quasars in general about a decade ago (see sec- 
tion |Ij). There has been great progress in reverberation mapping of Seyfert 1 
galaxies and low-luminosity quasars but nothing is known about line variability 
in high-luminosity quasars. Do high-luminosity quasars emission lines respond 
to continuum changes, as seen in Seyfert 1 galaxies and low-luminosity quasars? 
What is the relative amplitude of the response, if any? What is the lag of the 
response, reflecting the light-travel time across the BLR? Do high-luminosity 
quasars BLR sizes scale with AGNs luminosity, and lie on a continues relation 
from the faintest Seyferts to the brightest quasars? 

The best way to answer these questions is by monitoring spectrophotomet- 
rically a sample of high-luminosity quasars. However, the difficulties in carrying 
out such a program are enormous. As those objects have apparent magnitude 
much fainter than the PG quasars one needs a much larger telescope to monitor 
them (such as 8m class telescope). Also, if the continuum variations and the 
BLR size are scaling with the luminosity, much longer monitoring periods are 
required (of order of 10 years - though some preliminary results can be obtained 
after only few years of observations). 

In spite of the difficulties we initiated a program to monitor high-luminosity 
quasars. Our sample consist of 11 quasars which were chosen to have high 
northern declination, redshift of 2 to 3.4 (to include the Lya and Civ UV lines 
in the optical region), and observed V magnitudes in the range of 16-18 mag. 
The sample is being monitored photometrically in B and R bands each month 
at the WO since 1994 November. Preliminary results are presented in Figure H. 
All 11 quasar show variations in the R-band flux of ;^ 0.1 mag. The quasar 
S5 0836+71 has the largest variability which is ~ 0.3 mag. The variation of 
the high-luminosity quasars are smaller than the variations found for the PG 
quasars by about a factor of 5 in magnitude. While in a comparable time-period 
the variations in the PG quasars were in the range of 0.5-1 mag (Giveon et al. 
1999) the typical variations of the high-luminosity quasars are only 0.1 mag. 
However, while the monitoring period is comparable in the observer time frame 
the rest-frame monitoring period of the high-luminosity quasars is about a factor 
of 4 smaller than the PG quasars monitoring period and only amounts to ~ 1.5 
years. 

Spectrophotometric observations are needed in order to check a correspond- 
ing variations in the emission lines. We carried out few preliminary observations 
toward the brightest objects in the sample using the WO and SO telescopes. 
Our best observations are demonstrated in Figure 0. We present two observa- 
tion epochs separated by two years, each with total integration time of 4 hours. 
The observing technique is, as described in section |2], of using a comparison star 
in the slit simultaneously with the quasar. During the 2 years period a contin- 
uum variation of about 10% is seen. This corresponds to the fading this object 
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Figure 5. R-band photometry light curves for 11 high-luminosity 
high-redshift quasars. All objects show variations of ;^ 0.1 mag. 
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Figure 6. Two spectra of HS 1946+7658 {z = 3.051, V = 15.8 mag) 
taken on 1997 July at SO 2.3 m telescope (highest spectrum) and on 
1999 July at WO Im telescope (dotted line spectrum). The difference 
spectrum is shown at the bottom. Total integration time is 4 hours 
for each spectrum. Though continuum variation is evident, no line 
variations were found. 

shown in the R-band photometry (Figure q). No line variations is detected in 
the difference of the two spectra. 

Figure |6| demonstrates that monitoring high-luminosity high-redshift quasars 
is feasible and interesting results can be obtained, however, it also demonstrates 
the difficulty to carry out such a monitoring using small to medium size tele- 
scopes. Even though the total integration time is 4 hours and the observing 
conditions are good the S/N of the spectra is poor. It is clear that a large 
telescope is needed to carry out this spectroscopic program. Early in 2000 we 
started monitoring a subsample of our 11 quasars with the 8 m Hobby-Eberly 
telescope which is partly owned by the Pennsylvania State University. We hope 
that within a short time we will be able to present preliminary results from this 
long term monitoring program. 



6. Summary 

We reviewed the final results from a spectrophotometric monitoring of a large, 
optically selected quasar sample (for the detailed study see Kaspi et al. 2000). 
We find time lags between the optical continuum and the Balmer-line light curves 
for all AGNs with adequate sampling. Our work doubled the number of AGNs 
with a measured time lag (i.e., BLR size). We also increased the available lu- 
minosity range for studying the size-mass-luminosity relations in AGNs from 
two to four orders of magnitude allowing, for the first time, to construct reliable 
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relations. Using all AGNs with known BLR size we find that the BLR size scales 
with the 5100 A luminosity as ^0-70±o.03^ This is significantly different from pre- 
vious studies and is in contradiction with simple theoretical expectations, both 
suggesting i?BLR cc L^'^. We also obtained a mass-luminosity relation for AGNs, 
M oc L^-^^^-^ ^ which, however, has a large intrinsic scatter. These findings im- 
pose new and strict limitations on theoretical considerations and constrain any 
theoretical models which attempt to explain the AGN phenomenon. 

One of our major current goals is to further increase our knowledge about 
continuum and line variations and BLR size in high-luminosity quasars (L > 

10 ergs s^^). To that end we are initiating a program to monitor a sample of 

11 high-luminosity high-redshift quasars. We presented preliminary results of 
the continuum variations in those objects to be ^ 0.1 mag. Spectrophotomet- 
ric observations of the sample is on the way and we hope that with 8 m class 
telescopes such as the Hobby-Eberly telescope we will be able in a few years to 
extend our knowledge about BLR size to the most luminous quasars. The con- 
clusion of this long-term monitoring will hopefully complete the reverberation 
studies to cover the entire AGNs luminosity range of 10^^-10^^ ergss"^. 
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